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Abstract. We construct one soliton solutions for the nonlinear Schrodinger equation with variable 
quadratic Hamiltonians in a unified form by taking advantage of a complete (super) integrability 
of generalized harmonic oscillators. The soliton wave evolution in external fields with variable qua- 
dratic potentials is totally determined by the linear problem, like motion of a classical particle with 
' acceleration, and the (self-similar) soliton shape is due to a subtle balance between the linear Hamil- 

, tonian (dispersion and potential) and nonlinearity in the Schrodinger equation by the standards of 

£SJ ' soliton theory. Most linear (hypergeometric, Bessel) and a few nonlinear (Jacobian elliptic, sec- 

ond Painleve transcendental) classical special functions of mathematical physics are linked together 
' (-h '' through these solutions, thus providing a variety of nonlinear integrable cases. Examples include 

O |. bright and dark solitons, and Jacobi elliptic and second Painleve transcendental solutions for several 

L variable Hamiltonians that are important for current research in nonlinear optics and Bose-Einstein 

^3 , condensation. The Feshbach resonance matter wave soliton management is briefly discussed from 

this new perspective. 

> • 

"sj - ! 1. Introduction 

O . 
IT) ■ 

f\) ■ Advances of the past decades in nonlinear optics, Bose-Einstein condensates, propagation of 
soliton waves in plasma physics and in other fields of nonlinear science have involved a detailed 
study of nonlinear Schrodinger equations (see, for example, [9J, [17] > [SI], |149j . |151] and references 
therein). In the theory of Bose-Einstein condensation [12], |114j . from a general point of view, the 
dynamics of gases of cooled atoms in a magnetic trap at very low temperatures can be described by 
an effective equation for the condensate wave function known as the Gross-Pit aevskii (or nonlinear 
Schrodinger) equation [16], [69], [ID], [ID], |108] and [113] . Experimental observations of dark 
and bright solitons [15] , [20] . [H], [7_5] and bright soliton trains [9], |129j . |130] in the presence of 
harmonic confinement have generated considerable research interest in this area [IB] . [58J. 

The propagation of an optical pulse in a real fiber is also well described by a nonlinear Schrodinger 
equation for the envelope of wave functions travelling inside the fiber [1], [7J, [T7J, [ID], [SI], [77]. A 
class of self-similar solutions that exists for physically realistic dispersion and nonlinearity profiles 
in a fiber with anomalous group velocity dispersion is discussed in [IS], [ID], [97] , [98], [US], |123] . 
[125] . which suggests, among other things, a method of pulse compression and a model of steady- 
state asynchronous laser mode locking [98]. Solutions of a nonhomogeneous Schrodinger equation 
are also known for propagation of soliton waves in plasma physics [13], [22], [23], |104j . 
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Integration techniques of the nonlinear Schrodinger equation include Painleve analysis (10], [17] . 
[3D], ED], 132, [S3], E3], [SI], [HE], US], Hirota method [63], R [SI], Lax method [TU], [SI], 
[El], [151] . Miura transformation [93] . [95], inverse scattering transform and Hamiltonian approach 
0, [3], [6], [55], [59], [TDH] among others [21], [IS], [SI], [S3], pJJ7], [ED]. Although the classical 
soliton concept was developed for nonlinear autonomous dispersive systems with time being an 
independent variable only, not appearing in the nonlinear evolution equations (see [126] , [127] for 
highlighting this point), connections between autonomous and nonautonomous Schrodinger equa- 
tions have been discussed in [2], [3J, [27], [62], [82], [98], [112] . [116] and [155] (see Remark 2 for an 
explicit transformation). The formation of matter wave solitons in Bose-Einstein condensation by 
magnetically tuning the interatomic interaction near the Feshbach resonance provides an example 
of nonautonomous systems that are currently under investigation [T6], [SS], [130J . 



We elaborate on results of recent papers [9], [TU], [TJ, [T7], [52], [SS], [23], [62], [SS], [IS], [79], 
S3, U HIS], [123], [121], [US], [126], [127], [129], [HQ], [136], [U7J, [US], [119] on construction 




of exact solutions of the nonlinear Schrodinger equation with variable quadratic Hamiltonians (see 
also [136] and [151] , [152] . [153] ). In this paper, a unified form of these soliton-like (self-similar) 
solutions is presented, thus combining progress of the soliton theory with a complete integrability 
of generalized harmonic oscillators. We show, in general, that the soliton evolution in external 
fields described by variable quadratic potentials is totally determined by the linear problem, similar 
to the motion of a classical particle with acceleration, while the original soliton shape is due to a 
delicate balance between the linear Hamiltonian (dispersion and potential) and nonlinearity in the 
Schrodinger equation according to basic principles of the soliton theory. Examples include bright 
and dark solitons, and Jacobi elliptic and Painleve II transcendental solutions for solitary wave 
profiles, which are important in nonlinear optics [1], [22], [23], [7_S], [79], [97J, [123] . [125] . [151] and 
Bose-Einstein condensation [12], [9], [121], [129], [138] . 

The paper is organized as follows. We present a unified form of one soliton solutions with 
integrability conditions, and sketch the proof in the next two sections, respectively. In Section 4, 
more details are provided and some simple examples are discussed. Section 5 deals with a Feshbach 
resonance management of matter wave solitons. In the last section, an extension of our method is 
given and a classical example of accelerating soliton in a linearly inhomogeneous plasma [22], [23] is 
revisited from a new perspective. An attempt to collect most relevant bibliography is made but in 
view of a rich history □ and the very high publication rate in these research areas we must apologize 
in advance if some important papers are missing. 



where the variable Hamiltonian if is a quadratic form of operators p = —id/dx and x, namely, 




2. Soliton-Like Solutions 




Hip + gip + h\ip\ 2 ip, 



(2.1) 



iipt = —a (t) if> xx + b (t) x 2 ip — ic (t) xif) x — id(t)if> + g (x, t) if) + h (t) \ip\ 2 if> 



(2.2) 



Ref. [45] presents a detailed source on classical papers in the soliton theory. 
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(a, b, c, d are suitable real-valued functions of time only) has the following soliton-like solutions 

V>0M) = ^===ex P (i(a(t)x 2 + fi(t)xy + 1 (t)y 2 )) (2.3) 
xF(/3 (t)x + 2 1 (t)y) 

((p is a real constant, y is a parameter and /x, a, /3, 7 are real- valued functions of time only given 
by equations (I2.16l) - (l2.22p below), provided that 

g = g a (t) fi 2 (t) (fi (t) x + 2 1 (t) y) m , h = h a (t) fi 2 (t) /i (t) (2.4) 

(go and ho are constants and m = 0, 1). As we shall see in the next section, these (integrability) 
conditions control the balance between the linear Hamiltonian (dispersion and potential) and non- 
linearity in the Schrodinger equation (12. 2p thus making possible an existence of the soliton-like 
solution (with damping or amplification) in the presence of variable quadratic potentials (see also 
[12] . [63] . |126j . |127] and |154] for discussion of important special cases; Remark 2 provides an im- 
portant interpretation of relations (12 .4p as a complete integrability condition for the nonautonomous 
nonlinear Schrodinger equation (12. ip when m = 0). 

Here, the soliton profile function F (z) of a single travelling wave-type argument z = fix + 2 r yy 
satisfies the ordinary nonlinear differential equation of the form 

F"(z)=g z m F(z) + h F 3 (z). (2.5) 

If m = 0, with the help of an integrating factor, 

— ) =C + g F 2 + -h F 4 (C is a constant) , (2.6) 
dz J 2 

which can be solved in terms of Jacobian elliptic functions [8], [53], [81], [143] . When m = 1, 
equation (12. 5 j) leads to Painleve II transcendents [6], [28], [31], [33], [8T] . 

The variable phase is given in terms of solutions of the following system of ordinary differential 
equations: 

^ + b + 2ca + Aaa 2 = 0, (2.7) 
at 

dB 

-j- + (c + Aaa) fi = 0, (2.8) 

-1 + a fi 2 = (2.9) 
dt 

(see Ref. [31] and the next section for more details), where the standard substitution 

a * (2.10) 
4a (t) n (t) 2a (t) v 1 

reduces the Riccati equation (12 .7p to the second order linear equation 

At" - r (t) n' + 4tr (t) n = (2.11) 

with 

a! , / n , , „ d f a! d'\ , . 

r(t) = 2c + 4(i, a(t) =a6-cd + d 2 + - . (2.12) 

a 2 \a d ) 

(Relations with the corresponding Ehrenfest theorem for the linear Hamiltonian are discussed in 

Ref. [36].) 
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It is worth noting that in the soliton-like solution under consideration (12. 3p linear and nonlinear 
factors are essentially separated, namely, the nonlinear part is represented only by the profile func- 
tion F of a single travelling wave variable z = 8x + 2^y as solution of the nonlinear equation (12. 51) . 
Letting 0x + 2^y =constant, one obtains 

8' 

x' + ^x = 2a8y (2.13) 



a' 



and 

for the soliton velocity and acceleration with the aid of (12. 9j) . Then, by (I2.7f) — ( T278T) : 

x" - -x' + (lab - c 2 + c (- - - H x = (2.15) 
a \ \a c J J 

that is similar to equation of motion of a classical particle (damped parametric oscillations). 

The initial value problem for the system (I2.7p - (l2.9p . which corresponds to the linear Schrodinger 
equation with a variable quadratic Hamiltonian (generalized harmonic oscillators [15], [19], |60j, 
[145] . |15Uj ). can be explicitly solved in terms of solutions of our characteristic equation (12. lip as 
follows pi], [3B], [132], [133]: 

H (t) = 2fi (0) tM> (t) (a (0) + 70 (*)) , (2.16) 

ait ^ ao{t) -WWT^y (2 - 17) 

PK) 2(a(0)+ 7o (*)) W 1 ^ 

•*> = ^- 4(a(KU) - (2 ' 19) 



where 



a ° {t) = Mt) Mt) ~ Mt) ' (2 ' 20) 

AW = -^|, A(t) = exp(- f\c(s)-2d(s)) ds), (2.21) 



n, m 1 M*) , ^(0) 

7o(^) = ^ 77TT 7TT + TTyT ( 2 - 22 ) 

2yUi (0) /x (t) 2a (0) 

provided that /i and /ii are the standard solutions of equation (12. lip corresponding to the following 
initial conditions /x (0) = 0, fi' (0) = 2a (0) ^ and y, x (0) ^ 0, (0) = 0. (Formulas tfZTM - fl2~22l) 
correspond to Green's function of generalized harmonic oscillators; see, for example, [31], [36], [50] . 
[88] , [132] , [133] and references therein for more details.) 

The continuity with respect to initial data, 

lim a(t)=a (0) , lim 8{t)=8 (0) , lim 7 (*) = 7 (0) , (2.23) 

t— >0+ t— >-0+ t— >0+ 
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has been established in |132] for suitable smooth coefficients of the linear Schrodinger equation. 
Thus the soliton-like solution (I2.3P evolves to the future t > starting from the following initial 

data: 

ip(x,0) = limVOM) (2-24) 
t->o+ 



e i<t> 



; exp (i (a (0) x 2 + /3 (0) xy + 7 (0) y 2 )) 



xF(/3 (0)x + 2 7 (0)y), 

where 0, /1 (0) , a (0) , /3 (0) , 7 (0) and y are arbitrary real parameters (see also (16.211) for a more 
general solution of this form). 

Remark 1. When m = 0, the gauge transformation ip = \ { x , t) exp [ig (7 (t) — 7 (0))] changes the 
original equation Ii2.2\) into 

iXt = -a (t) Xxx + b (t) x 2 x - ic (t) x\ x -id(t)x + h (t) \x\ 2 X, (2-25) 
where a, b, c, d are suitable real-valued functions of time t only and 

h(t) = h of 3 2 (0)f, 2 (0) a{t) f, {t \ (2.26) 

which is more common in practice. Once again, classical solution of the linear equation Ii2.11\) . 
namely, our characteristic function fi (t) , completely controls the specific form of the nonlinearity 
factor h (t) required for creation of the soliton (an extension is given in Section 6; see also Refs. [12] . 
[126] and |127j for important special cases). 

Remark 2. A simple change of variables, 

1) = ± = e u ** X (Px,i), (2-27) 

transforms the nonautonomous equation ^2.2\) with conditions \2.J$ , when m = 0, into a standard 
autonomous nonlinear Schrodinger equation with respect to the new variables £ = fix and r = 7 : 

iXt + goX + h \x\ 2 X = X(t, (2-28) 

which is completely integrable by advanced methods of the soliton theory [2], [6], [ST], |151j (see also 
[45] and references cited in the introduction). This observation provides an alternative approach 
to derivation of our equations \2. 7| )-( fIT70j) . An extension of the transformation \2.21ty is given in 

m- ' 

3. Sketch of the Proof 

Following [34] (see also [22], [HI] and [118] ). we are looking for exact solutions of the form 

^ = A(x,t)e iS ^ t \ S{x,t) = a{t)x 2 + P{t)xy + j{t)y 2 (3.1) 

(y is a parameter). Substituting into (12.21) and taking the imaginary part, 

A t + ((4aa + c) x + 2a(3y) A x + (2aa + d) A = 0. (3.2) 

For the real part, equating coefficients of all admissible powers of x m y n with m + n = 2, one gets 
our system of ordinary differential equations (I2.7I) - (I2.9I) of the corresponding linear Schrodinger 



6 



ERWIN SUAZO AND SERGEI K. SUSLOV 



equation with the unique solution (I2.16p -f l2.22p already obtained in Refs. [31], |132j . |133j and/or 
elsewhere. In addition, an auxiliary nonlinear equation of the form 

aA xx = gA + hA 3 (3.3) 

appears as a contribution from the last two terms. With the help of (12. 8 p and (I2.10p our equation 
(13. 2 p can be rewritten as 

A t - (^x - 2a^y\ A x + )^-A = 0. (3.4) 
Looking for a travelling wave solution with damping or amplification: 

A = A(x,t) = — f== F (z) , z = co (t) x + a (t) y, (3.5) 

y/fj, (t) 

one gets 

c' x + c[y = Ojx - 2af3y^j c (3.6) 
with Co = and C\ = 27 (or z = (3x + 2jy). Then equation (13. 3p takes the form 

T~2 F ( z ) = ~% F (*) + -^r F " W . ( 3 - 7 ) 

az z ap z ap z jj, 

which must have all coefficients depending on z only in order to preserve a self-similar profile of the 
travelling wave with damping or amplification. This results in the required equation (12. 5p under 
the balancing conditions (12 .4p and our proof is complete. (An extension is given in Section 6.) 

Remark 3. Assuming in ft3.4\ ) that 

A(x,t) = -±=B(£, T ), (3.8) 

y/fj, (t) 

where £ = /3x and r = 2jy, implies = B T with a general solution B (£, r) = F (^ + r) = 
F((3x + 2 iy ). 

4. Details and Examples 

A brief description of the method under consideration is as follows. In order to obtain soliton-like 
solutions (12. 3p explicitly, say in terms of elementary and/or transcendental functions, one has to 
solve, in general, the nonlinear equation (12.51) for the profile function F (z) in terms of Jacobian 
elliptic functions [5], [53], [SI], |115j . |143] (some elementary solutions are also available), when 
m = 0, or in terms of Painleve II transcendents, when m — 1 (it is known that if m > 1, this 
equation does not have the Painleve property [6], [81]). In addition, one has to solve the linear 
characteristic equation (12.111) . which has a variety of solutions in terms of elementary and special 
(hypergeometric, Bessel) functions [IT] . [87] . [103j . |115j . [142] . Many elementary solutions of the 
corresponding linear Schrodinger equation for generalized harmonic oscillators are known explicitly 
(see, for example, [33],[SS],[2S],[37],[5U],[SB], [133], [H5], [150] and references therein). Then, the 
linear part allows determination of the travelling wave argument z = (3x + 2 / ~yy and the damping 
(or amplifying) factor /i" 1 / 2 of the soliton-like solution (12.31) . Our balancing conditions (12 .4p control 
dispersion, potential and nonlinearity in the original nonlinear Schrodinger equation (12.21) . which is 
crucial for the soliton existence. (An extension is discussed in Section 6.) 
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4.1. Nonlinear Part. When m — 0, equation ( 12. 5 p is integrated to the first order equation (12. 6p 
and (the corresponding initial value problem) can be solved by the reduction of elliptic integrals 
in terms of Jacobian or Weierstrass (doubly) periodic elliptic functions [E], [53], [79], [81]. We are 
interested in real-valued solutions. Some of the classical nonlinear wave configurations are given by 

F(z) = | 9o + V9 2 o-^C h ^ ' (41) 



x 



cn (#o - 2C /i ) 



1/4 



9o + V9o - 2C Q h 
2y/gf-2C h 



l/2> 



if hn < and 



F(z) = I -9o + V(ti-2C <M (4 2) 



/?o 



xsn 



Coho \ ( go — \f g 2 ~~ 2Cq/io 



l/2> 



if go < 0. Here, cn(w, k) and sn(w, /c) are the Jacobi elliptic functions [8], [53], |143j . Familiar special 
cases include the bright soliton: 



-/i cosh (y^o^j 



with Co = in ( 14. ip and the dark soliton: 



fW= ^ t anh(y^| (4.4) 

with Co = (?o/ (2^o ) i n ( 14. 2p . when cn(u, 1) = 1/ coshu and sn(w, 1) = tanhu, respectively (the real 
period tends to infinity). More details can be found in Refs. [8], [53], [57J, [79], [143] . |151] and/or 
elsewhere. 

If m = 1, the substitution F (z) = g)j 3 ^j2~JTu i w (() and £ = -z^c/ 3 transforms (12. 5p into the 
second Painleve equation, 

w" = (w + 2w 3 . (4.5) 

In the limit w — > 0, this equation reduces to the Airy equation and its solution may be thought 
of as a nonlinear generalization of an Airy function [5]. There is a one-parameter family of real 
solutions w = Ak (C) that are bounded for all real ( with the following asymptotic properties: 

kAi (C) , C -> +oo 
A k(0 = { ^ a. ( 4 -6) 



r (CP 174 sin (s (C) - 0„) + o (|Cr 1/4 ) , C -oo 



Here, Ai(£) is the Airy function, —1 < k < 1 provided 7^ 0, r 2 = — 7r 1 In (1 — k 2 ) , 

^(C) = ^IC| 3/2 -^ 2 ln|C| (4.7) 
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and 



6q = —r In 2 + arg T 



1 r 



+ 4* 1 



2sign (fc)) 



(4.8) 



These asymptotics were found in [5], |122j and, eventually, had been proven rigorously in 
(see [6], [14], [28], [29], [31], [33], |135] and references therein for study of this nonlinear Airy function; 
graphs of these functions are presented in [28], see Figure 1 for an example: w,^ = Ai/ 2 (x) .) 
An application to a soliton moving with a constant velocity in linearly inhomogeneous plasma is 
discussed in Section 6. 




** X 



Figure 1. Graphs of the Airy function .5Ai(x) and the nonlinear Airy function 
u>. 5 = A1/2 (x) (red and blue in an electronic version, respectively, [2"8]). 



4.2. Linear Part. Generalized harmonic oscillators [15], [49], [60], |145] , P-50J, which correspond 
to the Schrodinger equation with variable quadratic Hamiltonians, are well studied in quantum 
mechanics (see also [34], [35], [50], [83], [88], |133] and references therein for a general approach 
and known elementary and transcendental solutions). A few examples include the Caldirola-Kanai 
Hamiltonian of the quantum damped oscillator [T9] . |43j . [72], |141] and some of its natural mod- 
ifications; a modified oscillator considered by Meiler, Cordero-Soto and Suslov [92], [2Z], and the 
degenerate parametric oscillator |38j ; the quantum damped oscillator of Chrusciriski and Jurkowski 
[26] . and a quantum-modified parametric oscillator among others. Green's functions are derived in 
a united way in Ref. |36j . 

4.3. Examples. Combination of linear and nonlinear parts by our formula f)2.3p results in numerous 
explicit soliton-like solutions for corresponding nonlinear Schrodinger equations. It is worth noting 
that in this approach most linear and some nonlinear classical special functions of mathematical 
physics are linked together through these solutions. 
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4.3.1. Nonlinear Optics. In the simplest form, 



one gets [3 



a U) = — ^ — p (t) = — ^ — (4.io) 

W l-4a * l-4a t 



/3 2 t 

7 (*) = 7o + : - A (i (t) = /io (1 - 4a t) 

1 — Aa t 



(/io, ao> A), 7o are constants) and 

(3 x + 2 (tq + (/3 2 - 4a 7o) t) y 
1 - 4a t 



(4.11) 



g(x,t) = - »i\, 2 z m (m = 0,l), (4.12) 
(1 - 4a t) 

»W = -^- (4-13) 

(Traditionally, a = and m = with ^ = x ex P (iQoPot) [81], [151] . The case m = 1 is discussed 
in Section 6.) 



The case b = c = 0, 



iipt = -aip xx - idip + gtp + h \ip\ 2 ip, (4.14) 



is of interest in fiber optics (see, for example, @], [7], [61], [77], [78], [79], [97], [98], [ITS], [TTS], p3] - 
[124] . [126j . |127] and references therein). Here, all parameters a(t) , d(t) and ft, (£) are functions 
of the propagation distance t = z and this equation describes the amplification or attenuation (if 
d is positive) of pulses propagating nonlinearly in a single-mode optical fiber, where ip (t, x) is the 
complex envelope of the electrical field in a comoving frame, x = r is the retarded time, a (t) is the 
group velocity dispersion parameter, d (t) is the dispersion gain or loss function, and h (£) is the 
nonlinearity parameter [75], [79]. 

The substitution if; = x e_A > A (t) = J ' d (s) ds results in 

■ d * , , t -2A | |2 /. 1K x 

l ~dt = ~ a dx^ + 9X+ x ' ^ ^ 

which, of course, can be solved by the method under consideration, but a standard change of the 
time variable, 

»t 

a(s) ds, (4.16) 

o 

transforms this equation into the previous one. Just replace t — > r (see also [78], [79], |123j . |124] and 
[125] , where this simple observation has been omitted). More general transformations are discussed 
in Refs. [2], [27J, [35], [82], [98], [UJ, [IIS], [133], [155] . 
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4.3.2. Harmonic Solitons. In a similar fashion, one can show that the nonlinear Schrodinger equa- 
tion of the form 

* = o + ^ X + 777 777^ X X (4.17) 

at 2 \ ox 2 J 2 (cos £ + 2a sin £) 

has the following explicit solution: 

X (x,t) = —= F(z), (4.18) 

V |/io| (cost + 2a sin£) 



where 

z _ Pox + (270 cos t - (fig - 4a 7 ) sin £) 3/ 
cos £ + 2a sin £ 

and 

2a cos £ — sin £ 2 Ai^y 



(4.19) 



S(a;,£) = 7 U + ^ (4.20) 

2 (cos £ + 2ao sin £) cos £ + 2a sin £ 

270 cos £- (/3^ -4a 7o)sin£ 2 go/? 2 sin £ 



y + 



2 (cos £ + 2a sin t) 2 (cos£ + 2a sin£) 

provided 

F" = g F + h F 3 (4.21) 

(m = 0, /i 7^ 0, a , /3 , 70, <7o> ^0 an d y are arbitrary real constants). The reader may choose the 
profile function F, say in one of the forms (I4.ip - fl4.4p . 

If m — 1, the Schrodinger equation, 

M W d 2 ^ 2 A , ffo/3 2 , . 

2-7T- = - — TT^7 + X W M (4.22) 

9£ 2 V dx 2 V 2(cos£ + 2« sin£) 2 ^ V 7 



2 (cos £ + 2ao sin £) 

has a solution-like solution of the form (I4.18P (last term in f)4.20p should be omitted) provided that 

F" = g zF + h F 3 } (4.23) 

which solution is given in terms of the nonlinear Airy function (() . Further details are left to 
the reader (see also Section 6). 



5. Matter Wave Solitons 

5.1. Gross— Pitaevskii equation. Discovery of Bose-Einstein condensates in ultra-cold gases of 
weakly interacting alkali-metal atoms has stimulated intensive studies of nonlinear matter waves on 
a macroscopic scale (see, for example, [TB], [2H], [SB], [13], [13U] ). The Gross-Pitaevskii equation for 
a zero-temperature condensate of atoms, confined in a cylindrical trap V (x, y) = muj\ [x 2 + y 2 ) /2, 
and a time-dependent harmonic confinement, which can be either attractive or expulsive, along the 
z direction, V t (z,t) = mJ 2 (£) z 2 /2, is given by [12], [22], [M], [85], P33, M, [T27] : 

= (-|1a + VU (r, t) + U | V (r,tf + iM>) * (r, i) , (5.1) 
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where U = ATrfi 2 a s /m, a s is the s-wave scattering length, m is the mass of the atom, V ext = 
Vq (x, y) + V\ (z, t) and the condensate interaction with the normal atomic cloud through three- 
body interaction is phenomenologically incorporated by a gain or loss term r\ (t) . If the interaction 
energy of atoms is much less that the kinetic energy in the transverse direction, then the substitution 



, 1 / . x 2 + y 2 A (t) \ 

* (r ' t)= v^ exp r^-^f + 2 1 (5 ' 2) 



xip ( — , u±t ) , A (t) = / rj (s) ds 

\a± J Jo 



allows one to reduce the three-dimensional Gross-Pitaevskii equation ( 15. ip to the following one- 
dimensional nonlinear Schrodinger equation in new dimensionless units ( = z/a± and r = oj±t : 

^ = I + ^ (?) CV) + k (r) H 2 ip. (5.3) 



dr 2 V d( 

Here, 



K(r) = 2e A ^, c 2 (r) = 4, a± = \f^~ (5-4) 

and ao is the Bohr radius (see Refs. [68], [76], [93], [99], |10Qj . [1 llj . |121] for more details). 
Letting 

"0 = Xexp (i / #(s) J (5.5) 



with the function g (r) given by (15. 9p below, equation ( 15. 3 p can be transformed into ( 12. 2p . where 

a — -, b = -u 2 (r) , c = d = 0, h = n{r). (5.6) 

Then our characteristic equation (12. lip take the form 

//" + w 2 (r) // = 0, (5.7) 

which describes the motion of a classical oscillator with variable frequency [BZj. (This equation 
coincides also with the Ehrenfest theorem for the corresponding linear Hamiltonian [36J.) Choosing 
the standard solutions /x (t) and /ii (r) with /x (0) = 0, p,' (0) = 1 and /ii (0) 7^ 0, f/ t (0) = 0, one 
can use formulas (I2.16l) - (l2.22p with c = d = in order to solve the linear problem in quadratures. 
This gives the soliton travelling wave variable z = (3x + 2 r yy and the following balancing conditions: 

/3 2 (0)/i 2 (0) 
2^ (r) 

2 (r) = f 9 (s) ds = g (7 (0) - 7 (r)) (5.9) 



when m = 0. 
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5.2. Feshbach Resonance. The properties of Bose-Einstein condensed gases can be strongly al- 
tered by tuning the external magnetic field. A Feshbach resonance management for Bose-Einstein 
condensates has been discussed from experimental and theoretical perspectives by many authors 
(see, for example, [I], [25], [41], pED], [39], [56], [63], [66], [67], [71], [74], [EH], [90], [HE], [T09] . 
[TTO] . [TT4] . [TT7] . [T27] . [T28] . [T30] . [T3T] . [138], [139], [154] and references therein). The Feshbach 
resonance is a scattering resonance in which pairs of free atoms are tuned via Zeeman effect into 
resonance with vibrational state of the diatomic molecule [39], |130j . |139j . (They are known as Fes- 
hbach resonance because of their similarity to scattering resonances described by Herman Feshbach 
in nuclear collisions.) The strength of the nonlinearity U is defined in terms of s-wave scattering 
length a s , namely, 

47r^ 2 a, 

U = -, (5.10) 

m 

and dependence of atomic collision cross section due to existence of the metastable state [32], |56j, 
[71] . |128] enables a s to be continuously tuned from positive to negative values. (The scattering 
length also determines the formation rate, the spectrum of collective excitations, the evolution of 
the condensate phase, the coupling with the noncondensed atoms, and other important properties 
|40j . |114j .) As follows from the experiments, the s-wave scattering length is the following function 
of the applied magnetic field [96] : 

a,(B) . (5 . n) 



a \ B — B 

(The Feshbach resonance provides, so to speak, a continuous knob to adjust the atom-atom inter- 
action from repulsive to attractive, and from weak to strong [130J. Thus it is possible to study 
strongly interacting, weakly or noninteracting, or collapsing condensates [71] . all with the same 
alkali species and experimental setup. When the nonlinearity U = 0, one deals with linear modes 
of a macroscopic harmonic oscillator [76]; see, for example, [S3], [50], [88J and references therein 
for a detailed treatment of the corresponding quantum oscillator with variable frequency.) In the 
empirically established expression (15. lip . -Bo is the resonant value of the magnetic field, o m is the 
off-resonance scattering length and parameter Ao represents the resonance width in units of the 
Bohr radius a (see [39], [96], |127j . |138j . |139] and references therein for more details). Feshbach 
resonances have been observed in 85 Rb at 164 G [41], [TT7J, [10], in 23 N a at 853 and 907 G [67j and 
have also been identified in 6 Li [66], [IDS], [129], [T30] . 



5.3. Matter Wave Soliton Management. The Feshbach resonance provides an effective practi- 
cal tool for experimental study of the matter wave solitons. Indeed, for creation of a certain soliton 
configuration one needs to satisfy the following condition: 

/3 2 (0)/i 2 (0) _ A / Ap \ 

h e = ^[1 + — -j (5.12) 

in order to synchronize the Feshbach resonance and harmonic trap. (Here, both sides have the 
same simple pole structure, which can be used in experimental setting.) This equation allows 
determination of the classical law of motion (kinematics in z direction) of the expectation value 
fi = r^) with respect to the linear part of Gross-Pitaevskii Hamiltonian ( 15. ip . when U = 77 = 0, 
in terms of a suitable applied magnetic field B near the Feshbach resonance. (The synchronized 
harmonic trap oscillation frequency should be found from the classical equation of motion (15. 7p as 



SOLITON-LIKE SOLUTIONS 



13 



B = B + A + — — ~. ; N 2 . (5.15) 



B = B + A + — — ^ | 2 , (5.16) 



a; 2 = — Vice versa, the required tuning magnetic field is given by 

R=R 4a oo A e A /x 

^° + 4 aoo A e V - M 2 (0) /i 2 (0) ' { ° 6> 

if a particular law of motion /i is obtain by integration (dynamics) of the classical equation (15. 7p . 
Our criteria of the wave matter soliton management are consistent with ones obtained in Refs. [12] , 
|63j . |126j . |127J and |154j . if the classical equation of motion (15 ,7p is taken into account (this point 
seems not emphasized in these papers). 

5.4. Examples. Harmonic matter wave solitons, which correspond to uo 2 = uoq/uo 2 ^ =constant in 
the nonlinear Schrodinger equation (15. 3p . namely, 

f = \(~W + ^) +2 >^ ^° (5 ' 14) 

can be produced in Bose-Einstein condensates by tuning the external magnetic field near the 
Feshbach resonance as follows 

uh O fi o f3 2 A 
4aoo (2«o sin out + us cos cot) — uhofioP 2 , 

Letting uj — )■ 0, one gets 

h fi /3 2 A 
Acioo (2a T + 1) - h Q n Q Pl 

when ojq = (a similar case has been recently discussed in Refs. [63] and |127j ; it is of interest to 
analyze possible experimental setup; see also |154j ). 

The reader may find more details on the synchronization of Feshbach resonance and harmonic 
trap, explicit soliton configurations, and available numerical and experimental results in recent 
papers [12], [63], [126], [127], [148] (see also [9], [10], [33], [52], [78], [79], [140], [149], [151] and 
references therein). 

6. Generalization 

If an arbitrary linear combination of operators p = —id/dx and x is added to the quadratic 
Hamiltonian in equation (12. 2p . namely, 

itpt = -a (t) ipxx + b (t) x 2 ip - ic (t) xip x -id(t)ip - f (t) xip + ig (t) ip x + h (t) \ip\ 2 tp, (6.1) 

one can look for exact solutions in a more general form 

i> = A{x,t)e iS{x ^\ (6.2) 

S (x, t) = a (t) x 2 + /3 (t) xy + 7 (t) y 2 + 5 (t) x + e (t)y + n(t) + £ (t) 

(y is a parameter, we are separating contributions from linear k and nonlinear £ parts in the constant 
term). The linear part has been already solved in [51] and |132j . One has additional equations 

7 r 

— + (c + 4aa)5 = f + 2ag, (6.3) 
^ = (g - 2aS) 0, (6.4) 
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d,K 



^ = gS- a5 2 (6.5) 
at 



to the system (I2.7p - fl2.9p . whose solutions are given by 



x m r m M*) (^(Q) + go (*)) , . 

Ht) = 6o(t) ~ 2(a(0) + 70 (t)) ' (6 ' 6) 

m /3(Q) (^(Q) + go (*)) , fi7 ^ 

£ (*) = £ — --— , 6.7 

2 (a (0) +7o(t)) 

«<«> - "(°> + *m- Sct - (6 - 8) 



Here, 



5o(t) 



A(t) 

/"o (t) y 



a (s) / 2a (s 

2a (t) A (t) A* a(s)a( S )A (a) 

Mo(*) io (// (s)) 



f/,S 

M5 



(6.9) 



eo(t) = -^7^^o(t) + 8/ ^iA^ (/ , o(s)5o(s)) ds (6 . 10 ) 



Jo M s ) V a(s) J 

(*) = (t) - 4 f 00 * ( S )) 2 d* (6.11) 

-2 (/.o (.) 6 (*)) (f (s) - i&g (s)) ds 

Jo M s ) \ a \ s ) J 

with 5 (0) = -e„ (0) =g(0)/ (2a (0)) and k (0) = (see Refs. [34] and p2] for more details). 
Our equation (13. 2 p takes the form 

At + ((4aa + c) x + 2a/3y + 25a -g)A x + (2aa + d) A = (6.12) 
and the nonlinear equation (13. 3 j) becomes 

aA xx = ^A + /iA 3 . (6.13) 
at 

Once again, with the help of (I2.8p . (I2.10p and (16 .4p . equation (I6.12p can be rewritten as 

A t - (jx - 2a(5y + f)A x + ^A = (6.14) 

and looking for a travelling wave solution of the form 

A — A(x,t) — —j== F (z) , z = c {t)x + c 1 (t)y + c 2 (t), (6.15) 
(*) 

one gets 

c' x + c[y + c' 2 = I — x-2a(3y + -\ c (6.16) 
with c = (3, C\ = 27 and c 2 = e (or z = (3x + 2^y + e). Then 

£f W = ^ W + (6,7, 
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and our balancing conditions are given by 

§ = g a (t) d 2 (t), h = h a (t) d 2 (t) p (t) = h o 2 (0) /i 2 (0) a ®*® (6.18) 
at /i (t) 

with £ = g (7 (0) — 7 (£)) according to ( 12.9ft . For the soliton velocity, 

3' 

x' + jx = 2a/3y + 2ab - g, (6.19) 
thus extending ( 12.131) . Equation of motion is given by 

x" - ^-x' + ^Aab -J + c^-^yjx (6.20) 

= 2af+(^-c\g- g' 

as a nonhomogeneous generalization of (12. 151) (forced damped parametric oscillators). 
As a final result, our solitary wave solution has the form 

ip(x,t) = exp (i (ax 2 + (3xy + 7 (y 2 — go) + 6x + ey + k) ) (6-21) 

xF(0x + 2-fy + e), 

where the elliptic function F satisfies the nonlinear equation (12.51) with m = and 0, y, go and ho are 
real parameters. Time- dependent functions fi(t) , a (t) , B (t) , 7 (t) , 6 (t) and n (t) are given by our 
equations (I2.16l) -( l2.22l) and ( 16.61) — (16. Ill) (as in the corresponding solution of the linear Schrodinger 
equation [31], [132] ). 

Example 1. A soliton motion with acceleration in linearly inhomogeneous plasma was discovered 
in Refs. [22] and [23] (see also [13], [137] ). For a modified equation, 

where k, h , a , (3 an d ^0 are constants, we get // (t) = /i (1 + 4a t) and 

a (t) = — ^ — , B (t) = — ^ , (6.23) 

l + 4a t V ; l + 4a t 

Pot r/^ , So + kt 



1 + 4a t 1 + 4a t 

28ot(5 + kt) k 2 t 3 t(5 + kt) 2 

£ ( t ) = £ o 1 _i_ /I -a ' Kt = K o-^ r— - 

1 + 4a r 3 1 + 4a t 



with 



z = Bx + 2jy + e (6.24) 
R x-2t {Boy + 6 + kt) 
= TTi^t + 2 ^ + ^o 

in our particular solution ( 16 .21)) . The classical case [22], [23] corresponds to a = and ho^ofio = ~ 2 
(with A; — > —k). The reader may choose the profile F in one of the forms (I4.ip - (l4.4p . 
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Example 2. A similar case occur, if one takes m — 1 in (14.121) . The corresponding Schrodinger 
equation is given by 

M + ** = 9oft , + Wo 2 ^ (6 25) 



one gets 



and 



Here, 



where 

= /3 x + 2 ( 7o - (gg - 4qq7q) t) y 
1 + -In,,/ ' 

With the help of the gauge transformation, 

^-**xM, i- ^-«y , <«*> 

ofa.^ = «qZ 2 + /to + (To ~ (go ~ 4a 7o) t) V 2 (Q 3Q) 

yOMo 1 + 4a t y 

and the soliton profile is defined (as a solution of the second Painleve equation) in terms of the 
nonlinear Airy function (C) with asymptotics given by (14.60 — (I4.8|) . (Graph of one of these 
functions, = Ai/ 2 , is presented on Figure 1 from |28J.) It is worth noting that, in a contrast to 
the previous case, our A-soliton moves with a constant velocity when a® = 0. Further details are 
left to the reader. 
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